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1 Introduction 



The phenomenological study of the structure functions of deep inelastic scattering (DIS) 
is well developed [Q and its role in future colhders can hardly be overestimated 0]. 
One of the gaps in the state-of-the-art QCD description of the structure functions is the 
region of large Bjorken x, the so-called semi-inclusive region [|l], ^. The two limitations 
are insufficient data and lack of pure theoretical understanding of this region. The 
QCD analysis shows large perturbative and non-perturbative corrections which make 
the parton model inapplicable at large x. 

In general, the Operator Product Expansion (OPE) ||^ provides a well-established 
framework for analyzing the structure functions in QCD. It allows one to expand the 
moments of the structure functions Fa, where a = L, 2, in inverse powers of the hard 
scale l/Q"^, 

^ t=q,g ^ 

(1) 

expressing each term in the expansion as a product of the matrix elements of certain 
local composite twist— n operators between hadronic states, {Oj^\N)), and the corre- 
sponding coefficient functions Cl ^(A^, a^)- The latter can be systematically calculated in 
perturbative QCD, while the former bring a set of universal non-perturbative parameters, 
corresponding to generalized moments of some non-perturbative distributions. 

The number of the matrix elements entering the expression for the A^-th moment 
of the structure function depends on the twist. The leading, twist-two flavour non- 
singlet, contribution to (|l|) depends, for arbitrary A^, on just one non-perturbative pa- 
rameter, the A^-th moment of quark distribution within the hadron. The flavour singlet 
contribution depends on the A^-th moment of both the quark and the gluon distributions. 
In any case, twist-two depends strictly on single particle distributions. The asymptotic 
behaviour of the structure function at large x corresponds, in moment space, to large 
A^ ~ 1/(1 — x). An additional simplification occurs in this limit at twist two: the contri- 
bution of gluons to the structure functions becomes suppressed as compared with that 
of the quarks 

Going over to higher twists, one encounters a completely different situation: the 
number of matrix elements entering the twist-n {n > 4) contribution depends on the 
moment A^ and it grows rapidly with A^. The reason is that the high-twist contribution 
depends on novel multi-particle distributions, which describe the correlations between 



partons in the nucleon and which do not have a simple partonic interpretation @]-[^ 
For different moments A^, the high-twist matrix elements measure specific projections of 
the multi-parton distributions. One consequence is that the dependence of the high- 
twist contribution to the structure functions is not described by the DGLAP evolution 
equations. 

Examining the large-x limit one finds that the twist expansion breaks down. The 
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physical reason for this is that for large-x kinematics the final state in the current 
fragmentation region is bound to be a narrow jet with a square invariant mass W"^ = 
(5^(1 — x)/x such that W"^ <^ Q^. In such a situation, the long-distance effects of the jet 
formation dominate, in spite of the momentum transfer being large: the twist expan- 
sion goes effectively in inverse powers of the smallest scale, 1/W'^. As indicated by recent 
phenomenological analysis |jl2[-|jT5|], the problem becomes acute when W'^ ^ few GeV^. 

The breakdown of the OPE at large x is three-fold: the perturbative expansion breaks 
down, the power expansion breaks down, and the number of non-perturbative parameters 
increases. Perturbative corrections to the twist-two coefficient functions C2^{N,as) is 
enhanced at large by Sudakov logarithms 

Cla{N, a.) ~ E «r(«- 1^2- N + h^ In^™-! N + --). (2) 

m>0 

These large corrections originate from an incomplete cancellation between the contribu- 
tion of real and virtual soft gluons at large x. To obtain a reliable perturbative description 
of the coefficient functions, they need to be resummed to all orders in Ug- Thanks to 
the factorization property of soft radiation such resummation can be performed [rB|-[TU| 



and it takes the form of exponentiation. One expects to find similar Sudakov corrections 
in the coefficient functions of higher twist, C^„(A^, a^) with n > 2. Since higher-twist 
corrections scale as powers of l/W"^ ~ N/Q"^ rather than as terms of arbitrarily 

high twist in (|l|) become relevant at large A^. Consequently, the number of matrix ele- 
ments involved increases sharply. This means that, without additional insight into the 
properties of the high twists, the OPE parameterization of non-perturbative corrections 
to the structure functions in the semi-inclusive region becomes impractical and their 
scale dependence is obscure. 

In the absence of a practical OPE-based parametrization of power corrections to 
structure functions, various models were developed p0|-p9| and applied |0-||15|. In 



particular, the renormalon-based model became popular. Through the resummation 
of the factorially divergent perturbative series for the twist-two coefficient functions, 
infrared renormalons generate an ambiguous higher-twist contribution. The renormalon 
approach to power corrections is based on the consistency of the theory as well as that 



of the OPE: all ambiguities in the resummed perturbative calculation must eventually 
cancel out. Therefore, the presence of ambiguities ~ /i^^/Q^" in the twist-two coefficient 
functions implies the existence of genuine power corrections of the same form, that is of 
the twist— (2n-|-2). In this way infrared renormalons in the twist-two coefficient functions 
detect some part of the higher twist contribution. 

The main assumption made by the renormalon model (e.g. [^-[^) is that the 



power corrections accessible to renormalons dominate. A physical argument supporting 
this assumption was never provided. A further assumption which, in our view, makes 
this model inapplicable at large x, is that multi-gluon contributions are discarded by 



considering the large-/3o limit PDI: in this approximation one resums a set of diagrams 



corresponding to a single dressed gluon. Since at large x multiple gluon emission is crucial 
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at the perturbative level it is unlikely that power corrections would be entirely 
associated with a single emission. Therefore, in contrast with the above mentioned 
renormalon model, we shall not base the parametrization of power corrections at large x 
on the single dressed gluon. In our approach multiple gluon emission plays a major role 
also at the non-perturbative level. 

The first part of this paper is devoted to show how the renormalon ambiguity cancels 
out within the OPE, leading to unambiguous predictions to the structure functions. The 
conjectured cancellation of infrared renormalons in physical observables is fundamental 
to our understanding of power corrections in QCD (see e.g. Nevertheless, not much 

has been done to test it. In this context the example of deep inelastic scattering plays 
an important role: since the DIS structure functions admit an OPE, the mechanism 
of cancellation can be checked by explicit calculations. The case of the longitudinal 
structure function Fl{x) = Fi{x)—F2{x) was analyzed in this spirit in Ref. ||20|, It was 



shown that the renormalon ambiguity in the twist-two coefficient function indeed cancels 
against that of the quadratic ultraviolet divergence of the twist-four matrix elements. 
The cancellation of renormalon ambiguities for the F2{x) structure functionQby the same 
mechanism will be demonstrated here. 

The OPE is based on the separation of short and long-distance effects into the co- 
efficient functions and hadronic matrix elements, respectively, by introducing the fac- 
torization scale /i. This scale provides an infrared cut-off for momenta of particles in 
Cl^^{N,as) and, at the same time, it sets up the ultraviolet renormalization scale for 
{01"'\n)). Although both the coefficient functions and the matrix elements have a 
nontrivial dependence on /i, this dependence cancels in the expression for the structure 
functions (|l|). In particular, varying fi one finds that the twist-four operators mix among 
themselves at the level of logarithmic ultraviolet divergences, ~ In/x^, as well as with 
the operators of twist two at the level of quadratic ultraviolet divergences, ~ /x^. The 
former divergences give rise to DGLAP evolution equations, while the latter introduce 
an ambiguity in separating the contribution of twist two and twist four. This ambiguity 
is compensated in (|^) by ~ yU^/Q^ contribution to the twist-two coefficient functions in- 
duced by infrared renormalons. Thus, the cancellation of renormalon ambiguities follows 
from the independence of the DIS structure functions of the factorization scale. 

Let us now return to the renormalon model and address the meaning of the underlying 
assumption, that is the dominance of those power corrections which are associated with 
infrared renormalons. The predictive power of the model relies on the assertion that the 
dependence of the renormalon ambiguity on the Bjorken variable x, or equivalently on 
the moment A^, represents that of the full higher-twist contribution. This is realized if the 
contribution of the higher-twist matrix elements is dominated by the part that has the 
highest possible ultraviolet divergence, and thus mixes under renormalization with the 

^The ambiguity cancellation mechanism is believed to be general. In particular, it should apply also 
to the parity-violating structure function F3 as well as the spin structure functions. Note, however, 
that explicit expressions for the twist four contribution, equivalent to p^ , are not yet available in these 
cases. 



3 



leading twist. When viewed in the context of the OPE, the "renormalon dominance" 
assumption is therefore understood as |3^ "uhraviolet dominance". It was used 



in [|25|-[|28[, and [0^0] to develop the phenomenology of power correction to the 
DIS structure functions. 

To leading order in the large-/5o limit, i.e. at the level of a single dressed gluon, 
the ambiguities in the twist-two coefficient functions of both Fi{x) and ^2(0;) appear as 
power corrections ~ /^^/Q^ and ~ i^^lQ^ only. In the OPE, they are compensated by 
the contribution of the operators of twists four and six, respectively, that are given by 
the convolution of a specific function of x with the twist two parton distributions. Then, 
the "renormalon dominance" assumption implies that 

• Contributions to the structure functions of twist higher than six are suppressed, 
and can be ignored. 

• The dominant twist-four and twist-six contributions are proportional to the twist- 
two distribution. Therefore, these contributions can be parameterized by a single 
non-perturbative parameter at each twist which controls the overall normalization. 

Because of the small number of parameters, the renormalon-based model for power 
corrections seems to be very appealing phenomenologically. It is important, however, 
to keep in mind the strong assumptions involved. A priori, from the OPE point-of- 
view, ultraviolet dominance seems hard to justify. Higher- twist matrix elements measure 
multi-parton correlations within the hadronic target. They encode information on the 
detailed structure of the hadron, which is inaccessible at twist two. The renormalon 
model reduces this rich structure to an overall normalization of the twist-two matrix 
elements. 

The main shortcoming of the renormalon model, which is particularly relevant at 
large x, can be understood through perturbation theory. The renormalon ambiguity 
described above emerges from a specific set of diagrams in which a single dressed gluon 
is emitted. These diagrams dominate the large-order behaviour of the series. However, 
it is well understood that at large x multiple gluon emission is important: due to the 
phase-space constraint, i.e. the small invariant mass of the hadronic system, W"^ <^ Q^, 
the emission of soft and collinear gluons is enhanced. In perturbation theory this results 
in the appearance of the Sudakov logarithms in the coefficient functions (j^). These log- 
arithms exponentiate in the Mellin transform (moment) space owing to the factorization 
property of soft and collinear gluons [|16|-||19|. This implies that the large- asymptotic 
behaviour of the perturbative coefficients is controlled by multiple emission from the 
primary hard parton (s). An entirely different set of diagrams than the single dressed 
gluon. 

Analyzing the properties of the resummed perturbative expressions, one finds that 



multiple soft gluon emission also affects the power corrections The systematic 



way see this is by considering the large-order behaviour of the Sudakov exponent |36 



It turns out that the coefficients of sub-leading Sudakov logarithms are enhanced 
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factorially due to infrared renormalons. This implies that any Sudakov resummation 
with a fixed logarithmic accuracy is insufficient unless W"^ ^ A^. In the large x limit 
one is obliged to resum all the logarithms that are factorially enhanced. Here infrared 
renormalons hit again: the resummed exponent is defined only up to a power-suppressed 
ambiguity. The ambiguity must be compensated by genuine power corrections, which 
exponentiate together with resummed perturbative expansion. 

The emerging structure of power corrections becomes quite different from the renor- 
malon model: due to multiple emission all powers of l/W"^ ~ N/Q"^ are present 



dxx''-'F,{x,Q') = ClaiN,a,) (Of (iV)) 



1 + 



ra>l 



(3) 



where a = 2,L and Q'^/N is ffxed as — oo, and terms which are suppressed by 1/N 
are neg lected. Here, (Of (A^)) and Cl^{N, as) are the twist-two quark-operator ma- 
trix elements (the gluon contribution is negligible) and the corresponding Sudakov- and 
renormalon-resummed coefficient functions, respectively, and k^^^'^A^" are A^-independent 
non-perturbative scales related to hadronic matrix elements of twist-2(ri + 1) operators 
producing the dominant contribution as x ^ 1. 

Eventually, in order to describe the behaviour of the structure function at large x, 
one has to resum an infinite set of power corrections on the scale Q'^/N. In the shape 
function approach 0, ^ the infinite set of scales K^f^A^, with n = 1,2,..., is replaced 
by a single non-perturbative function. The series on the r.h.s. of (|^) defines the Laplace 
transform of the shape function (see Eq. (^) below), which is a new non-perturbative 
distribution which governs the shape of the structure function at large x. The shape 
function approach has already proven useful in several processes, including Drell-Yan 
production 



_I3|, 0, and event-shape distributions pl|-p9|. 
Independently of the renormalon analysis, we examine in this paper the structure and 
the physical origin of the dominant power corrections at large x within the framework of 
the OPE. At first sight, the increasing number of local matrix elements of higher-twist 
operators in the large- A^ limit seems to be an obstacle to any use of the OPE to study 
this limit. Nevertheless, analyzing the twist-four contribution at leading order using its 
representation in terms of non-local light-cone operators p[-|Tl| we find that, as at the 
leading twist, a significant simplification occurs in this limit. 

Physically, this simplification occurs owing to the dominance of particular configu- 
rations in the final and initial states. The dominating final state turns out to be the 
same as at twist two. The initial state differs, of course: the quark-antiquark light-cone 
operator of twist two is replaced by a quark-gluon-antiquark operator at twist four. 
However, in the dominating configuration the gluon has a small longitudinal momentum 
fraction. Consequently, the corresponding twist-four matrix elements essentially depend 
only on the light-cone separation between the quark and antiquark, just as the twist-two 
matrix element. 

We find that the relevant correlation measured by twist four at large x corresponds to 
very specific configurations, which resemble twist two. This leads us to finally to the our 
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main conjecture, that the simplification of higher twist at large x amounts to ultraviolet 
dominance, namely the dominance of these operator matrix elements that mix under 
renormalization with the leading twist. Note that, in contrast with the renormalon 
model, we do not assume here the dominance of single gluon contributions; multi-gluon 
contributions are relevant. The ultraviolet dominance assumption implies that the twist 
two matrix element factorizes out of the twist expansion at large x. Moreover, for the 
twist four contribution to be totaly free of factorization-scale dependence, the resummed 
coefficient function must have the same logarithmic dependence on the factorization scale 
as the twist two coefficient function. Thus both the parton distribution function and the 
resummed coefficient function factorize, as summarized in (^). 

This factorization implies, in particular, that the large limit has a remarkable 
feature: a subset of matrix elements which dominate the twist four contribution decou- 
ples from the rest and satisfies a closed set of evolution equations. In fact, such an 
approximation is known for the twist-three contribution to the spin structure functions 
of deep inelastic scattering and we plan to extend it to twist four in a forthcoming 
publication. 

We proceed as follows: in Section 2 we recall the standard formulation of the OPE, 
and then present the known results concerning renormalon ambiguities in the twist- 
two coefficient functions. In Section 3 we follow |]1T| describing twist four in terms 
of non-local light-cone operators and calculate their coefficient functions in moment 
space. We also examine the large-x limit and show that the asymptotic behaviour of the 
twist four contribution in this limit is associated with specific final states. In Section 
4 we verify explicitly the cancellation [ED], ^^ between ~ infrared renormalon 



ambiguities at twist two and the ambiguity in the definition of twist four matrix elements 
for both Fl and F2, to leading order in the flavour expansion. In Section 5 we return to 
the large x limit and show that the physical picture emerging from general considerations 
at twist four is consistent with the ultraviolet dominance assumption. Our conclusions 
are summarised in Section 6. 



2 Twist-two analysis and infrared renormalons 

The spin-averaged structure functions, Fl{x) and F^^x), of deep inelastic scattering of a 
virtual photon with momentum q from a hadron with momentum p are deflned in terms 
of the hadronic tensor 

with X = Q^/2pq and = —q^ > 0. The latter is expressed through the imaginary 
part of the Fourier transform of the time-ordered (T) product of two electromagnetic 
currents, symmetrized over the Lorentz indices: 

W,.{x,q') = ^ Im {/ d'ye''y{p\zT [j^^{y)j^y{-y)] (5) 
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with jn{x) = eg^I/g(a;)7^\l/g(x) and Nf being the number of quark flavours^ At 
large the dominant contribution to (j^) comes from the vicinity of the hght-cone 
~ — > 0. This allows one to expand the T-product into the series of local com- 
posite operators with increasing twist (= dimension — spin), and, therefore, decreasing 
singularity for small y^. A more efficient (although equivalent) way of writing the OPE 
is based on expanding the T-product in terms of non-local light-cone operators of in- 
creasing twist |Tl[|. These operators serve as generating functionals for the infinite tower 



of high-twist local operators. This formulation of the OPE proves to be more convenient 
for our purpose. 

To leading order in the coupling constant, the twist-two contribution to the OPE of 
the symmetric part of the T-product in (|^) is given for — >■ by [|TT|] 



iT 



-1 tw-2 1 



X / du 







l + \\nuy'dl + 0{y' 



^{uy)i^{-uy) , (6) 



where d = d/dy and the sum over Nf massless fiavours is assumed. A path-ordered 
exponential between —y and y is required to make this operator gauge invariant. We do 
not write it explicitly here. 

Going over to hadronic matrix elements, we introduce the standard twist-two quark 
distribution g(^) by the relation 

{p\^{y)0{-y)\p), = 2{py) J^e^'^^^y^i^, f^') , (7) 

where y^ = and /i defines the normalization scale of the non-local light-cone operator. 
Inserting (|^) into and calculating the hadronic tensor, we arrive at the well-known 
expressions for the Born approximation to the twist-two part of the structure functions 

Tf (x) = , Tf ) ix)=xY.el [q{x) + g(x)] , (8) 

where q{x) = q{—x) stands for the antiquark distribution. Going beyond the leading 
order in (|]), one finds that the coefficient function of the non-local quark operator 
is modified by 0{as) corrections and, moreover, the T-product receives contribution 
from non-local light-cone gluon operators. The general expression for the twist-two 
contribution to the structure functions is then given by 

-Tr-2(x,g2) = /'f cL(e,«.(Q'), + (9) 

X J X t 



where a = L,2 stand for Fl and F2, respectively, and C2^a{^,as{Q^),Q'^/ fi'^) are the 
corresponding quark coefficient functions. The ellipses stand for the antiquark and the 
gluon contributions. 

^In what follows we shall assume for simplicity that quarks of different flavour have identical electric 
charges e,. 
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For our purpose, the gluon contribution to (|^) can be safely ignored for the following 
reasons. Our task in Section 4 is to verify the cancellation of infrared renormalons at twist 
two and ultraviolet contribution to the twist-four matrix elements, to the leading order of 



the flavour expansion. However, as was argued in ||2^, renormalons appear in the gluon 
matrix elements only at the next-to-leading order in the flavour expansion. Secondly, 
when examining the asymptotics of the structure function F2{x) in (y) for x — 1, one 
finds [^, ^ 1^, that the gluon coefficient function is subleading as compared with that 
of the quark. For example, to order the F2 coefficient functions are 



Clf{^,as{Q'),l) = 6{l-0 + CF- 

n 



' ln(l-0 \ 3/ 1 - 



(10) 



c|,2(e,«.(Q'),i) = iV/ 

TT 



^e(l - 2^ + 2^') [ln(l - - HO] '1+ 4^(1 - 



and the Fl coefficient functions are 

q'f(e,«.(Q^),i) = c^^e 

Cl,d^,o^s{Q%l) = Nf^2ai-0, (11) 

' TT 



where the factorization scale was set as yU^ = Q^. In both cases there is an additional 
suppression factor of (1 — ^) for the gluon coefficient function compared with that of 
the quark, and therefore the gluon contribution can be ignored provided that the gluon 
density is not anomalously large compared with the quark density for x — 1. We shall 
assume that this is indeed the case. Phenomenological studies shows that the gluon 
density decreases faster than that of the quarks, consistently with this assumption. 

To see intuitively why the hierarchy between the quark and gluon coefficient function 
at large x holds, one has to compare the final states described by the two contributions. 
For the quark contribution to F2{x) the final state consists of the outgoing energetic 
quark jet of small invariant mass Q'^{1 — x) surrounded by soft gluons. For the gluon 
contribution, coming necessarily through the box diagram, the cloud of soft particles 
is bound to contain a soft quark. Invoking the standard power-counting arguments, 
one finds that the contribution of a soft quark is suppressed by a power of the energy 
~ (1 — x) compared with that of a soft gluon. 

In addition one finds that owing to the different quantum numbers, the Fl coefficient 
functions are suppressed at large x by an overall factor of {1 — x), compared with F2. 
Taking the moments of (^ 

Fa{N, = r F,{x, Q2) a;A^-2rfa; = c|,,(iV, a,) (12) 
Jo 
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with qN = q{x) ^dx and C2{N) = Jq C2{x) x^ ^dx, one concludes that 

Fl-2(iV)/F*"-2(iV) ~ C2,LiN)/C2,2iN) = Oil/N) , (13) 

The vahdity of the hierarchy (^) proves to be independent of the twist. In particular, 
as we shall show below, the same hierarchy ( pTsD holds for the twist-four contribution, 
and it becomes a general property of the structure functions in the semi-inclusive region. 

The twist-two coefficient functions C2 as((5^), Q^Z/x^) have a well-defined pertur- 
bative expansion. To any order in the coupling constant, their perturbative expressions 
contain terms involving logarithmic dependence on the factorization scale fi as well as the 
terms regular as yU 0. The former are controlled by the DGLAP evolution equations, 
while the latter become sensitive to the behaviour of the coupling constant at low scale 
and increase factorially at high orders in (renormalons) . Infrared renormalons make 
the perturbative series Borel non-summable and thus induce ambiguous contributions to 
the coefficient functions suppressed by powers of l/Q^. 

The logarithmic dependence of the twist-two parton distributions q{x,fi) on the fac- 
torization scale fi reflects the In /i renormalization scale dependence of the corresponding 
operators (^. Through renormalization, different operators mix with each other. Sim- 
ilarly, higher twist n (with n = 4, 6, . . .) distributions admit power-like dependence on 
the factorization scale, ~ /i^". This reflects the fact that the higher-twist operators 
mix under renormalization with the leading twist. As a consequence, the separation of 
higher-twist matrix elements from those of the leading twist depends on the factoriza- 
tion scale. Obviously, the structure functions Fa{x, Q^) should not depend on this scale. 
Within the OPE (P, the logarithmic /i dependence cancels out between the coefficient 
functions and the matrix elements of the same twist, whereas the power-like /i^" de- 
pendence cancels between the contributions of different twists. From this we draw two 
important conclusions. The first is that in order to verify the cancellation of infrared 
renormalon ambiguities in the twist-two coefficient functions, one has to calculate the 
most divergent ultraviolet contribution to the matrix elements of higher-twist operators 
- the quadratic divergence for twist four, the quartic divergence for twist six, etc. The 
second conclusion has considerable implications on phenomenological applications: in 
general, higher-twist contributions cannot be consistently included^ without resumming 
the perturbative expansion (which is, by itself, ambiguous). The ambiguity cancels only 
when resummation is performed and power corrections are included. 

Renormalon contributions to the twist-two coefficient functions, ~ a"^^/?o ^!, can 
be extracted from the (unphysical, but technically useful) large Nf limit, since /3o = 
(IIC^ — 2Nf)/12 is linear in Nf. The leading order in the flavour expansion corresponds 
to the lowest-order diagrams in flgure where the internal gluon line is dressed by any 
number of fermion loops. 

To regularize a divergent renormalon contribution, we shall apply the Borel regular- 

•^A similar conclusion was reached in [ p3| . 
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Figure 1: Diagrams that contribute to the twist-two coefficient function at leading order. 



ization [gQ]. It amounts to representing the large- A^/ dressed-gluon propagator by 

and identifying singularities of the resulting momentum integral in the w-plane located 
at real positive u. The main advantage of the Borel method is that it allows one to 
regularize simultaneously the renormalon contribution to the coefficient functions and 
power-like ultraviolet divergences to the higher-twist matrix elements. 

The renormalon contribution to the twist-two coefficient functions, C22 and C|^, has 
been calculated in the single dressed gluon approximation in Refs. (see also |^ 

and [^]). It generates power-suppressed corrections to the coefficient functions of the 
form ~ l/Q^ and ~ only and, therefore, produce ambiguous contribution to the 

structure functions of twist four and twist six, respectively. In particular, the twist-four 
renormalon contribution has the following form 



f 

Jo 



dxx^-^A 



2x 

Fl 
2x 



tw— 4 



tw— 4 



4N 



N + 2 



+ 



12 



1 + 2 + 



-8-iV 



(14) 



with il){x) = d\nr{x)/dx and 7 being the Euler constant. Here, the notation was 
introduced for the dimensionless constant 6, defined by the integral 6 = j^du/il — u), 
which represents the difference between different deformations of the original Borel- 
integration contour [0 < u < 00) in the vicinity of the pole at m = 1. 
The expressions (|H) simplify significantly at large N, 



dxx^-^A 



dxx^~^A 



.2x 
Fl 



2x 



Cp NA^6 



tw— 4 



tw— 4 



% 



Qn 



(15) 
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so that the renormalon contribution to Fl{x,Q'^) is suppressed at large by a power 
of as compared with -F2(x, Q^). We recall that a similar property holds for the twist- 
two contribution, Eq. ([T3|) . In addition, Eq. (|15D suggests that, in agreement with our 
expectations (|^), the effective parameter of the twist expansion in the limit x ^ 1 is 
given by N/Q"^ rather than l/Q"^. An expression similar to (|l^) is available for the 
twist-six renormalon ambiguity In the Borel plane, this ambiguity appears as an 
additional pole, at u = 2. At large A^, this contribution contains an additional factor of 
N/Q"^ with respect to twist four. Again, this holds for both F2{x,Q'^) and Fl{x,Q'^). 

Consistency of the OPE of the structure functions requires that renormalon ambigui- 
ties (0) will cancel in (|ID against the ultraviolet divergent contribution to the twist-four 
matrix elements. For the longitudinal structure function, Q^), this cancellation 

was demonstrated in |2T]]. Below we will conduct a similar analysis for case of the 
structure function F2{x,Q'^), which according to (0) is expected to have a more sin- 
gular behaviour for x ^ 1. The cancellation of the twist-six renormalon ambiguities 
with ultraviolet divergent contribution to the twist-six matrix elements requires a com- 
plete analysis of twist-six operators and their coefficient functions. It is therefore more 
involved and will not be discussed here. 



3 Twist-four analysis 



To demonstrate the cancellation of the renormalon ambiguities (p!4D, we will need the 
expression for the twist-four part of the structure functions. This relies on the twist-four 
contribution to the symmetric part of the time-ordered product (^) . To leading order in 
the coupling constant it is given by [I 



iT 



i{^.{y)3u}{-y) 



tw— 4 



1287r^ 



(16) 



with two terms T^^^u and '72, representing the separate projections of the longitudinal 
and the transverse components, respectively, and the ellipses denoting the operators 
containing total derivatives which therefore having vanishing forward matrix elements. 
The operators Ta^f^u (with a = 2, L) are defined for — > as 



L, 



du 



U 



1 + Inn) Qiiuy) + uhiud^dy{\u.y'^ Qiiuy) 



;i7) 



and 



du{ - Iny^ df,dy 



u \nuQi{uy) 



u 



Q2{uy) 



+^{y^l^u + yud^ - g^y{y ■ dj) 



-m(1 + 3\nu) Qiiuy) + - Q2{uy) 



u 



;i8) 



Here, Qi{y) and Q2{y) are two specific combinations of non-local twist-four light-cone 
operators, which will be given later on. 
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Inserting (|T^, (0) and ( jTB]) into the expression for the hadronic tensor, (^), we 
parameterize, in analogy with (^, the matrix elements of the operators Qiiy) and Q2{y) 
in terms of the twist-four distributions -Di(0 and D2{C,), respectively. 



{pmy)\p) 



(19) 



with y"^ = 0. The symmetry of the T-product, Eq. (0), under y —y implies that 
Qii~y) = Qiiy)) and, as a consequence, -Di(0 = A(— (with i = 1,2). The functions 
Di{C,) defined is this way describe the momentum fraction distribution inside the hadron 
and vanish outside the interval — 1 < ^ < 1. However, in contrast to the twist-two dis- 
tributions, they do not have probabilistic interpretation and, in general, are not positive 
definite. As we will show below, they can be expanded into a sum over multi-parton 
configurations with definite weights, given by the twist-four coefficient functions. 

Calculating the twist-four contribution to the hadronic tensor according to (|^), we 
perform a Fourier transform of the forward matrix element of (|T^) and take its imaginary 
part. Matching the resulting expression into the definition of the structure functions, 
Eq. (^, we obtain after some algebra the following expressions for the twist-four contri- 
butions to the structure functions 



2x 
F2(x) 



2x 



tw— 4 



tw— 4 



^/;f«i+..f)D.(x/o. 

1 r /■! (if 

^ x{D2ix)-D,{x))+ -j^\5 + 6\nOD,ix/0 



(20) 



The expressions obtained for the twist-four contribution to the structure functions 
have a striking similarity to the twist-two contribution given by (^). Indeed, one can 
rewrite ( pO| ) in the same form as (|^) by replacing twist-two parton densities by their 
twist-four counterparts -Dj(0 and defining the corresponding coefficient functions C4^a(^). 
This analogy proves to be very useful as it reveals the important difference between the 
two terms entering the second relation in (pOl) : the first term where the distribution 
amplitude is evaluated at ^ = x, and the second, which involves an integral over the 
partonic momentum fraction ^. In the case of twist two, these two terms are associated 
with the contribution to the coefficient functions due to virtual and real perturbative 
corrections to the partonic cross section, respectively. The former corresponds to a 
single-parton final state, C2,a(0 ~ ^(1 ~ 0^ and it is singular at a; = 1. In general, 
to high orders in the perturbation expansion, the singular behaviour of the coefficient 
functions for x — 1 is associated with particular partonic (short- distance) final states 
in which the recoiling quark creates a narrow jet. Generalization of this analysis to the 
twist-four contribution is straightforward. 



Following ||T0[, we can rewrite the twist-four contribution to the structure function 
as the sum over all possible final states. Then, the two terms in the above-mentioned 
expression for F2{x) correspond to two different short-distance final states: a single 
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recoiled quark state and a state in which this quark is accompanied by an additional 
gluon carrying the momentum fraction 1 — ^. We also notice that the coefficient functions 
entering the expressions for Fl{x) and the second term in the expression for ^2(0;) in 
(pO|) are not singular as ^ — 1. As in the case of twist two, this suggests that the 
emitted gluon cannot be soft and, as a consequence, the corresponding contribution to 
the structure functions is subleading at large x compared with the first term in the 
second relation in (pO]). 

Indeed, let us examine the twist-four contributions (pOD in the large x limit. It is 
convenient to consider the expressions (^) in moment space 



dxx 



N- 



-1 Fl{x) 



2x 



1 1 



dxx 



N- 



^1 W 



2x 



tw-4 



tw-4 



4g2 {N + 2) 
1 



(21) 



16g2 



D2{N + 1) - 



(iV + 2)2 



where Di{N) = d^^^'^ A(0- Assuming that the moments Di{N + l) and D2{N + 1) 
have a similar behaviour at large N, we conclude from (^Tj) that the twist-four contri- 
bution to Fl{x) is suppressed by a power of N compared with that of -^2(0;). We recall 
that the same property holds for the twist-two contribution to the structure functions, 
Eq. (p!3D, including the renormalon contribution ([T5|). 

We also find that for large x, the twist-four contribution to -^2(0;) (this is not so for 
Fl{x)) has a natural approximation. 



F2{x) 



2x 



tw-4 



x{D2{x)-Di{x)) 



(22) 



Remarkably enough, this is similar to the expression for the twist-two structure function 
to the lowest order in a^, Eq. (|]). As explained above, the reason for this is that in both 
cases the partonic final states which determine the large-x asymptotics of the coefficient 
functions look alike: they correspond to the propagation of a single recoiled quark. 
Calculating high order corrections to the twist-four coefficient functions, we expect that, 
as in the twist-two situation, the large-x asymptotics will correspond to the propagation 
in the final state of narrow jet initiated by the recoiled quark. The difference between 
twists resides, at large x, in the different initial states. For the twist two, the initial state 
consists of a single parton and the coefficient function is given by the squared modulus of 
the diagonal transition amplitude, or, equivalently, the partonic cross section. For twist 
four, there can be one or more partons in the initial state, so that the coefficient function 
is given by the sum of diagonal transition amplitudes as well as their interference. This 
implies that twist-four coefficient functions do not have an interpretation of the partonic 
cross section but rather of correlations between different multi-parton states within the 



hadron 10 



Let us examine the properties of twist-four distributions Di{x) in more detail. Ac- 
cording to Eq. (|T^, these distributions are defined through hadronic matrix elements of 
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non-local light-cone twist-four operators Qi{y). Their expressions have been worked out 
| TT| and are given byQ 



m m- 



Qi{y) 



dv 



+ 4 / dt{l- vt) 0^{v, t- y) + {v- t) 0(,{v, t- y) 



+ ^ -y) 



Q2{y) 



dv 



-1 



4i(l - v"^) Oj{v- y)+4: dt\{2-v + t- 2vt) O^^v, t; y) 



+ {y^-y). (23) 



+ {v-t) 0,{v, t- y) + or y) + Or{v, t; y) 



Here, the notations for the "canonical" basis of light-cone operators Oi are according 
to p[. They can be classified according to their partonic content and are defined as 
follows. The quark-antiquark-gluon operators are 



0;{v-y) = -ge^^^^^ij{y)y^Yl^GP'^{vy)^{-y) 
Or{v;y) = g^{y)iy,D^G'^''{vy)^{-y) , 



(24) 



where 7/^ = and G^^^ = with being the color matrices in the quark represen- 

tation. Other operators include a quark, an antiquark and two gluons. 



0,{v,t;y) = g'^{y)y^G''\vy)y^G^^{ty)^^{-y) 

Oe{v,t-y) = zg\,p,^{y)yf'y^G"^^{vy)ypG^''{ty)j^j,^{-y), 

and C-even combination (sym) of two quark-antiquark pairs. 



(25) 



Onv,t-y) = '-^ 



Oriv,t;y) 



X 

4 

X 



'^{y)t''0{vy) - ^{vy)e0{y) 
'^{ty)e0{-y) - ^{-y)t^0{ty) 

'^{y)t^^^5^{vy) +^{vy)t^^-f5^{y) 
^(t?/)t'^^75*(-2/) + ^(-y)t"^75^'(t2/) 



(26) 



Inserting ( p3D into ([T9|), we can relate the distributions Di{^) to the matrix elements of 
the basis twist-four operators Oa{v;y), which, in turn, can be parameterized by intro- 
ducing the set of twist-four multi-parton distributions (see figure H) , 

{p\0^{v;y)\p) = 2{py) J d^,d^,e'^P^^^^^^'-^^+^^^'+^^^ U^,, ^2) 

{pmv;y)\p) = 2{pyf f d^,d^, e'^^y^^^^^'~^^-^^^^'+^^^ ^1,^2) • (27) 
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Figure 2: Forward matrix elements of the non-local light-cone operators 0^j{v\ y). 



Expanding both sides of (27) in powers of y we find that, as in the case of twist two, the 
moments of the distributions Tj(^i,^2) are related to the matrix elements of local twist- 
four operators. The distribution functions corresponding to the remaining operators 
from the canonical basis are defined similarly. 

The distributions T3j(^i, ^2) are the analogue of the twist-two distribution g(^) intro- 
duced in (|^). They describe the correlations inside the hadron between a state containing 
a quark and a gluon carrying the momentum fractions ^1, ^2 — ^1, respectively, and a 
state containing quark carrying momentum fraction ^2- These distributions vanish out- 
side the "physical" region —1 < ^1, ^2, — ^1 < 1- Using the hermiticity properties of 
the operators (p4|), [0^{v]y)f = —Os^—v; —y) and [0'j{v]y)'f = Oj{—v] —y) we deduce 
from (^) that 

T3*(ei,6) = ^3(6,^1) , T;(ei,6) = T7(6,6) • (28) 

In addition, time-reversal invariance of the matrix elements in (^) implies thatQ 

T3(6,ei) = T3(ei,6) , T7(6,6) = mi, ^2) . (29) 

From (^) and (^) one deduces that the distributions Ta(^i,^2) are in fact real-valued 
symmetric functions [p]] . 

Taking into account (pT]), it becomes straightforward to express the distributions 
Di(^) defined in (19) in terms of the canonical twist-four distributions Tj(^j). To this 



end one substitutes (|27|) into (|^) and matches the result of v integration into the r.h.s. 
of (|T^. Expanding the both sides of the resulting relation in powers of (py), we obtain 
the following expression for the odd moments of the distributions 

/' DiiO = (iV) + f 7(iV) + . . . 

Jo 

/' d^e-' ^2(0 = ■ fs{N) + 2 ■ fr{N) + ..., (30) 
Jo 



"'These expressions were taken from Ref. and a few misprints were corrected. 
^We thank O. Teryaev for useful discussion on this point. 
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where the elhpses stand for the contribution of the remaining distributions from the 
canonical basis, which we demonstrate in the next Section to be irrelevant for the can- 
cellation of renormalon ambiguities to leading order in the flavour expansion. Here we 
have introduced the notation for the generalized moments: 

f„(iV) = 161 c/«2Ta(ei,e2) $a,7V-2(ei,6), (31) 

with a = 3, 7. $3,Ar(^i) and ^7^N{C,i) are homogenous polynomials of degree N in the 
momentum fractions ^i. Inserting the expressions (PB| ) for the non-local twist-four oper- 
ators in terms of the canonical operators Oa into (p!9|) and (|30| ) uniquely determines the 
form of the polynomials ^a,N'- 



*3,iv(6,6) 



'^'7,iv(6,6) 



N+1 



S2 



N 



1 



k=0 



N-1 



N + 2 



%%'^3,iV+l(6,6) = 



k=0 



{k + l){N-k) 
N + 2 



kcN-k~l 



(32) 



Since the functions $3,Ar(^i) and ^'7,Ar(^i) are symmetric under interchange of and ^2, 



the moments (31) only receive nonzero contribution from the real part of the distributions 
Ti(^i), in virtue of (||). 

Repeating a similar calculation, one can obtain the expressions for the coefficient 
functions $Q,Ar(^j) (with a = 1, 2, 5, 6), defining the contribution of the remaining twist- 
four canonical operators to the r.h.s. of (pOl). 

Using (|21]) and ^ we get 



dx X 



N- 



-1 Fl{x) 



dxx 



N- 



2x 



-1 F2{x) 



2x 



tw— 4 



tw— 4 



X 



1 + 1 
4g2 {N + 2)2 

1 1 

16g2 {N + 2)2 



N{N - 1) T^{N + 1) + {N + 1){N + 8) Tj{N + 1) + 



(33) 



These are the general expressions for the contribution of the operators 0^{v\y) and 
0-j{v]y) to the moments of the structure functions. The moments Ta{N), Eq. (pl|), are 
given by the sum over matrix elements of local twist-four operators generated by the 
non-local light-cone operators, Oa{v;y). Each local operator enters into the sum with 
the weight given by the coefficients of the polynomials $a,Ar(^i), expanded in powers of 
momentum fractions. 



4 Cancellation of ambiguities 

The generalized moments entering the r.h.s. of (^) depend on the factorization scale 
yU that sets up an ultraviolet cut-off for the matrix elements of the twist-four opera- 
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tors Oa{v',y)- We now proceed to calculate the quadratically divergent, ~ /i^, contri- 
bution to the distributions Ta{C,i). According to (0), this induces an ambiguity in the 
definition of the twist-four structure functions, which should cancel against the infrared 
renormalon contribution, Eq. (|T5|). 

To match the renormalon calculation of twist two, we must calculate the ultraviolet 
contribution to the twist-four distributions in the same approximation, that is to leading 
order in the fiavour expansion. Since this contribution does not depend on the particular 
choice of the matrix element, it suffices to consider the matrix elements of the operators 
(j) \Oaiv;y)\p) between quark states, with the quark being off-shell, < 0, and y"^ = 0. 
To calculate the quadratic divergences we have to introduce a regularization. To be 
consistent with the twist-two calculation of the coefficient functions, we are obliged to 
use the same regularization, namely the Borel representation of the large-A^/ dressed- 
gluon propagator. The Borel variable u regularizes the ultraviolet divergent integrals 
over the gluon momenta and quadratic divergence amounts, in this regularization, to a 
pole at M = 1, leading to an ambiguity once the Borel integral is evaluated. 

Examining the matrix elements of the non-local light-cone operators, {p \Oa{v;y)\p), 
defined in equations (|2^)-(^6D, we find that, to leading order in the fiavour expansion, 
only the quark-gluon-antiquark operators 03{v;y) and Oi{v;y) (fig. ^ contribute to 
the quadratic divergence, and, as a consequence, mix with the twist-two operators. The 
relevant diagrams are shown in figures |^ and ^ 




Figure 3: Contributions to the quadratic ultraviolet divergence of twist-four matrix ele- 
ments. Insertion of fermion loops into the gluon lines is assumed. The contribution (b), 
corresponding to non-local operators with a quark-antiquark pair and two gluons, van- 
ishes. 

The irrelevance of the operators containing two gluons, 0^{y,t]y) and O^iy^t^y), 
was pointed out in for the case of Fl{x). It turns out that this holds also for ^2(0;). 



Let us briefiy explain why 05^Q{v,t]y) are not relevant. These operators mix with the 
leading twist only if both gluons are contracted. We recall that the twist-two renormalon 
contribution was calculated to leading order in the fiavour expansion, so we must apply 
the same approximation here. If both gluons are contracted to the quark or antiquark 
lines the contribution is sub leading in Nf and therefore irrelevant. There remains the 
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possible contribution where the gluons are contracted one to the other, as in figure ^ (b). 
However, as noted in pO[, this diagram vanishes for = 0, since there is no external 
scale. 

The vanishing of the contribution of the four-quark operators 0^2^ {v, t; y) of Eq. (^) 
is less obvious. The four quark operators (see figure || (a)) have the form 



Ci(si, S2; S3, S4) 

C^2(Sl, S2', S3, S4) 



^(Sl)t'^^(s2) ^(S3)t''^^(s4 



^(S3)t'^75^(s4) 



(34) 



In order for these operators to mix with the twist two operator (^, one of the quarks 
must be contracted to one of the antiquarks. In addition, for the contribution to be 
relevant in the large Nf limit a dressed gluon must be included. Contracting between a 
quark and an antiquark of different pairs, e.g. "^{32) and \E'(s3), and adding a dressed 
gluon with n fermion loops, the contribution is 0{g^{g'^N f)'"-). Being sub-leading in Nf, 
this contribution is irrelevant. On the other hand, contracting between a quark and 
an antiquark within a pair, e.g. ^^(si) and \E'(s2) as shown in Fig. ^ (b), an additional 
factor of Nf is gained so the contribution is 0{g'^Nf{g'^Nf)'^), which is relevant in the 
large Nf limit. Calculation of this diagram shows that {p |Oi,2(si, S2; S3, S4)|p) is indeed 




(a) 



(b) 



Figure 4: (a) The non-local four-quarks operators Oi^2', (b) a diagram contributing to 
the quadratic divergence of the matrix element of Oi,2(si, S2; S3, S4) in the large- A^/ limit. 



quadratically divergent. Nevertheless, because of the symmetry property 

(|}|Ci,2(Sl,S2;S3,S4)|p) = 6'i,2(p|Cl,2(s2,Si;S3,S4)b) (35) 

with $1 = —62 = 1, the symmetrized matrix element {p\OY^^{v,t;y)\p) is free of any 
quadratic divergence. As a consequence the four-quark operators do not contribute to 
the quadratic divergence in (p3D. 

To understand this property, it is useful to return to the construction of the OPE 
of eqs. (p^)-(|r^) as an expansion of local Feynman diagrams in powers of The 
four-quark operators Oif™(f , t; y) emerge from the diagram in figure |^ (a). The possibly 
relevant contributions of the non-local four-quark operators of figure ^ (b) correspond. 
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(a) 



(b) 



Figure 5: (a) Local Feynman diagram from which four-quark operators emerge by taking 
the operator product expansion; (b) The origin of figure § (b). 



before expansion, to the diagram in figure |^ (b). The latter vanishes since it contains a 
quark loop attached to one photon and two gluons, which is odd under (electric) charge 
conjugation. The symmetrized form of ( ^61) which nullifies its contribution is thus a 
manifestation of charge-conjugation invariance. 

We conclude that to leading order in the flavour expansion, the structure func- 
tions (^) receive a quadratically divergent contribution only from the generalized mo- 
ments T^j^N + 1), or equivalently from the distributions T3 7(^i,,^2)- The quadratically 
divergent contribution to these distributions, A^yT^ 7(^1,^2), was calculated in [^, 

Evaluating the diagram shown in Fig. |^ (a) with a dressed gluon and an off-shell 
external quark < at = 0, where the Feynman rules for the non-local operators 
are defined by (0), one gets 

AuvT3(6,6) = -^^''^ + 1)9(^1)^(^1 - 6) + (6 6)1 ^(6)^(6) 
A,vr7(ei,6) = s^^'*^ 11^(6)^(6 - 6) + (6 6)| ^(6)^(6) , (36) 

with the parameter of Borel regularization S defined in ( p^ ) and g(^) being the twist- 
two quark distribution function. Next, the ultraviolet ambiguity of T^jlN + 1) can 
be evaluated by substituting (|36|) into ( pi]) and performing the C,i integrals with the 
coefficient functions $3 7(^4) defined in (|32D . The result is 



AuvT3(iV+l) 
A^yfriN + 1) 



2V^(Ar + 1) + 27 - 1 + 



1 



N+1 



(37) 



N -2ip{N+ 1) -27 + 2 



+ 1 



Finally, substituting (^7^ into (|33| ) we calculate the ultraviolet ambiguity of the struc- 
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ture functions 



.2x 
Fl 



tw-4 



2x 



4/3o 



Ail){N + 1) +47 + 



+ 



12 



N 



Qn 



-4 4:poN + 2^'' Q 



(38) 



As expected, this ambiguity cancels exactly with the corresponding infrared renormalon 
ambiguity in the twist-two coefficient functions in Eq. (|T^. We see that when the 
twist-two coefficient functions are resummed and the (ultraviolet divergent part of) the 
twist four is included there remains no ambiguity at the level of 1/Q^ corrections in the 
OPE dll). 

We conclude that to leading order in the flavour expansion, infrared renormalons 
"probe" the contribution to the longitudinal, Fl{x,Q'^), and transverse structure func- 
tions, F2{x, Q^), coming from only two twist-four operators, O3 and O7, defined in (124 
An additional simplification occurs for x —>■ 1: from (^) it follows that for ^ 1 



A^yTriN + l)/A^yTsiN + 1) ~ AT, 



(39) 



and, therefore, the contribution of the operator O3 to the structure functions is subdom- 
inant at large x, as far as the cancellation of the renormalon ambiguity is concerned. 
As we will argue in the next Section, (^) is not limited the ultraviolet divergent contri- 
bution of the twist-four operators, but is a general property of these operators. Taking 
into account this hierarchy, we find that to leading order in the flavour expansion, the 
twist-four contribution to the structure functions for A^ ^ 00 is entirely associated with 
the quark-antiquark-gluon operator Oj, 



' N-i Fdx) 



2x 



1 



1 ^_iF2(x) 
axx 



2x 



tw-4 



tw-4 



16g2 



Tj{N + 1) + 0(A^°) 



Tj{N + 1) + 0{N^) 



(40) 



In addition, based on ( p^ ) we know that Tj{N + 1) contains a piece of the form 

Po 



(41) 



where we replaced the ambiguous 5 by a new A^-independent dimensionless parameter, 

[71 

Ki = 0{1), and, as before, we suppressed terms that are subleading at large A^. Under 
the ultraviolet dominance assumption the remaining piece is negligible, so 



T7(Ar + l) :^Tr(A^ + l) = K^^Nqj^A''. 

Po 



(42) 



This way the normalization of the twist-four contribution (^Dj) is fixed by a single non 
perturbative parameter k^^\ as anticipated in (Rf). 
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5 Large-x behaviour of twist four 



We saw that the ultraviolet dominance assumption leads to a great simplification of the 
parameterization of the twist-four contribution at large x: instead of a large number 
of local matrix elements, which increases N, this contribution is parametrized in terms 
of a single parameter. But is the ultraviolet dominance assumption justified? A clear 
cut answer to this question cannot be given within the framework of perturbative QCD, 
since it requires calculation of the matrix elements. However, by identifying the partonic 
configurations which become dominant at large x, we can understand the physical origin 
of the simplification that occurs at large x. 

We already saw, independently of the ultraviolet dominance assumption, that at 
large x, the twist-four contribution to F2 has a natural approximation (^2]) in terms of 
the distributions Di{C,), evaluated at ^ = x, which is analogous to the Born-level twist-two 
result At both twist two and twist four the recoiling system corresponds to a single 
quark. We note that for Fl at twist four, just as at twist two [|^, the recoiling system 
includes a hard quark and a hard gluon. In each case, at higher orders in perturbation 
theory the recoiling system is accompanied by soft and collinear radiation. The fact 
that the "bare" recoiling system is the same for twist four and for twist two opens up 
the possibility that the jet function is the same, as we shall argue below. Clearly, this 
condition is not sufficient. In general, the coefficient functions at twist four differ from 
those of twist two, since the latter strictly consist of single-particle initial states in the 
amplitude and the complex-conjugate amplitude, whereas the former contain more than 
one parton, at least in one of them. Indeed, -Dj(^), contrary to g(^), describes correlations 
of several partons in the hadron rather than a single-parton probability distribution. 

To investigate the relevant initial states probed at large x, let us apply ( P^D and 
to express the structure functions through the twist-four distributions Ta(6,6) defined 
in (pTf). In this way, we obtain the following expression for the longitudinal structure 
function in the momentum fraction representation. 



2x 



tw— 4 



4 

— / d^id^2 
J-i 



-mi, ^2) 



X 



X 



l + ln-le(6 



x) 



(43) 



+^7(6, 6) 



X 



(6 



6 - 6 + (36 - 6) In 



^(6 



x] 



(6 ^ 6 



and a similar expression for the transverse structure function. 



2x 



tw-4 



+T7(6,6 



x(5(6 



d^id^2 



x] 



T3(6,6 
x'ei^2 



X6{^2 



XI 



XI 



(6 - 6)^ (6 - 6)^el 



6-6 



36 - 56 + 6(3^ 




Here, the following notation was introduced for real-valued distribution functions, 



T3(6,6) 



2 [^3(6,6) 



T3(-6,-6 
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7^(^1,6) = \ [^7(6, 6) + T7(-ei, -6)] • (45) 



The coefficient functions entering (^31) and (|^) have poles at = ^2- Their regulariza- 
tion is uniquely fixed by going from the moment space in (|33|) to the momentum-fraction 
representation as 

1 - m,M^, (46) 



and it amounts to creating a punctur^ on the integration axis at = ^2- We recall 
that = ,^1 and S,q = —^2 have the meaning of the longitudinal momentum fractions 
carried by quark and antiquark, respectively, whereas = C,2 — C,i defines the momentum 
fraction of a gluon, so that + + = 0. The singularities of the integrand in (^) 
and (^^ at ^1 = ^2 can be attributed to a wee gluon in the initial state with a vanishing 
momentum fraction. As a check, it is straightforward to verify that the moments of (^) 
and (Q) coincide with (^), given the regularization (^). 

Obviously, Eqs. (|43| ) and ( PD have to be in agreement with ( ^OD as x — * 1. To see 
this, we note first that, owing to the presence of 6{C,i — x) and 9{^i — x) terms in the r.h.s. 
of (1^) and (^^, the momentum fraction of one of the quarks is restricted to x < < 1. 
Performing integration over the momentum fraction of the other quark, we notice that the 
dominant contribution comes from the vicinity of the pole at ^1 = ^2- As a consequence, 
the asymptotic behaviour of the structure functions ( |13| ) and ( ^^ is driven by the terms 
most singular at — ^2! 0. For Fl{x) these terms are ~ 6{^2 — ~ ■^2)^ and for 

F2{x) they look like ~ 5{^2 — — ^2)^ and ~ 6{^2 — a;)/(^i — ^2)^- Assuming that 

the non-perturbative distributions 73(^1,^2) and 77(^1,^2) have the same behaviour for 
^g = ^2 — ^\ ^ 0, we find that the large x asymptotics of (|43|) and (|^) is governed by 
the distribution amplitude T7(^i,^2)- 

The asymptotic behaviour of the structure functions can be associated with the prop- 
erties of the specific configuration probed at large x: it describes the correlations between 
a state consisting of a quark carrying momentum fraction ^1 ~ x and a state consisting 
of a quark and a gluon, carrying momentum fractions ^2 ^ x and = ^2 — G ^ 0, 
respectively. 

One can come to the same conclusion in the moment space by examining the proper- 
ties of the coefficient functions entering (PT|). Assuming that the functions Ta(^i, ^2) have 
a similar, non-singular behaviour, the large behaviour of Ta{N) is controlled by the 
large A^ behaviour of $o,Ar(^i, ^2)- One finds that at large A^ the functions $3,Ar((^i, ^2) 
and $7,Ar(^i, ^2), defined in the region —1 < ^1 , ^2 , ^1 — ^2 < 1, are peaked at the corners 
of the phase space ^1 = ^2 = =tl and, therefore, the moments of the structure functions 
F2{N) and Fl{N) receive dominant contribution at large A^ from the vicinity of these 
points. Comparing the contribution of the two operators in (pOD, we note that the peak 
value of $7,7v(Ci, 6) at ^1 = ^2 = ±1 is $7 = A^(A^^ - l)/6/(A^ + 2), which is larger than 

^The limits £ — > and iV — s- cx) do not commute. One has to calculate the moments first and then 
take the limit e ^ 0. 
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the corresponding peak of $3,7v = + 1. At large A^, the resulting hierarchy between 
the two is proportional to A^, in agreement with the conclusion we reached examining 
the ultraviolet divergent contribution in Eq. (|39D . 

We conclude that the configuration which dominates the Born-level twist four con- 
tribution at large x includes a gluon with a small longitudinal momentum fraction 

= ^2 — ii ^ ^- In this way it approaches the twist-two situation where the gluon is 
absent altogether. We also saw that the domination of the twist-four operator O-j with 
respect to O-^ at large x can be established independently of the ultraviolet dominance 
assumption. 

Finally, let us discuss the reason for the operator O-j to be dominant at large x as 
compared with O^. The difference between the two operators manifests itself in 



the different power of the pre-factor {py) in the r.h.s. of (0). The factor {py) has the 
meaning of the light-cone time. For x away from the end-point region x = 1 one has 
{py) = ^(1)- On the other hand, for x — ^ 1, or equivalently A^ oo, one finds that 
{py) ^ N . As a consequence, the contribution of the operator O-j to the moments of the 
structure functions is enhanced by a factor A^ as compared with O^. To show this, it 
proves convenient to return to (|2l|) and rewrite the moments Di{N + 1) at large (even) 
A^ in the following form 

A(A^ + 1)= / ^^a^A(0^ / rfee-^(^-«)A(0, (47) 
Jo Jo 

where we have replaced = exp(A^ln^) = exp(— A^ [1 — ^ + 0{{1 —C,)"^)]), and extended 
the integration region by taking into account that -Di(0 = for ^ > 1. Using (P^, we 
express -Dj(0 cis the Fourier transform of the matrix elements of non-local light-cone 
operators Qi{y), 

A(0 = r ^e-2^««(p|Q.(y)b) . (48) 

J-oo TT 



Substituting this relation into (^) and performing integration over ^ one obtains. 



To evaluate the integral over {py), one uses the fact {p\Qi{y)\p) has no singularities in 
the lower half planeQ. Then, assuming an appropriate behaviour of the matrix element 
{p\Qi{y)\p) infinity and closing the integration contour in (HSf) in the lower half plane 



we get the following remarkable relation [0] 



D,{N + l) = -{pU{y)\p)e 



-2i(py) 



(50) 

2{py)=-iN 



which is valid at large A^. Here the factor exp(— 2z(p?/)) compensates the phases coming 
from the wave functions of quarks separated along on the light-cone at the large distance 



^This can be seen from examining Eq. ( pl| ) for ^ > 1, where must vanish, so that closing the 

{py) contour in the lower half plane encloses no poles. 
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~ (py). Substituting ( |50[ ) into (|2T|), we find that the moments of the twist-four structure 
functions -Fl(x) and F2{x) at large N are determined by the matrix elements of non- 
local hght-cone operators {p\Qi{y)\p) analytically continued into the complex plane and 
evaluated at at large light-cone separations, 2{py) = —iN. A similar relation can be 
established for the large N asymptotics of twist two, 



(51) 

2{py)=-iN 



Recall now that the operators Qi{y) are given by a linear combination (p3D of twist- 
four light-cone operators, Oa- Comparing the two matrix elements defined in (0), we 
observe that {p\Ot{v; y)\p) is enhanced by an additional power of (py) as compared with 
{pIOs^v; y)\p). According to (|5DD, the contribution of the latter operator to the moments 
of the twist-four structure functions Fl{x) and ^2(0;) is suppressed by an additional 
factor with respect to the former. To arrive at (^) one has to take into account he 
fact that the v integrals^ over the position of gluon field on the light-cone are dominated 
hj V l/{py) ~ This brings the additional factor in front of {p\Oa\p) in 

agreement with large N behaviour of infrared renormalon contribution (W^). 



6 Conclusions and discussion 

6.1 The simplification of twist four at large x as ultraviolet 
dominance 

The breakdown of the twist expansion at large x, and, in particular, the increase in the 
number of local matrix elements required to parametrize the contribution of each twist, 
makes that framework, as it stands, impractical. Nevertheless, the twist expansion can 
be used, as we did here at twist four, to identify the source of the dominant contributions 
at large x, which scale as ~ Nh? /Q"^. 

On very general grounds, one expects some simplification in the description of DIS 
structure functions at large x. The leading twist simplifies in this limit since one is 
essentially probing only the valence parton distribution in the hadron. On the other hand, 
in this limit there are large perturbative and non-perturbative corrections associated with 
the formation of a narrow jet in the final state. Since the formation of the jet proceeds 
through colour exchange with the remnants of the target, the same corrections can be 
associated, by the OPE, to specific hadronic matrix elements. The physical picture we 
obtained is that the valence parton distribution remains the relevant physical quantity at 
large x, in spite of the fact that multi-parton Fock states give important contributions. 



^The role of the v integration differs for the case where \py\ is large from the expansion in small (py) 
used to derive Eq. (132). It should be noted that the function of v multiplying O7 in Eq. (ESh is related 



to the additional power of {py) in Eq. (27) by the requirement that the matrix elements, {p\Qi{y)\p), 
should be analytic at small {py). 
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which lead to the breakdown of the twist expansion. The multi-parton correlations 
measured by individual higher twist terms reduce in this limit to a trivial part which 
is associated with the power-like evolution of the valence quark distribution. Since thi s 
evolution is a consequence of the exchange of soft gluons (with momenta ~ ^Q'^{1 — x)) 
the finer structure of the hadron is not resolved. 

Our main conjecture in this paper is that the simplification of structure functions at 
large x is realised within the twist expansion through ultraviolet dominance, namely that 
each twist is dominated by the contribution which mixes under renormalization with the 
leading twist. To prove this statement, it is necessary to calculate the relevant matrix 
elements, a task which is definitely beyond reach of the tools we are using. Nevertheless, 
examining in some detail the twist-four contribution at leading order provides significant 
evidence supporting this conjecture. 

After reviewing the standard twist-two and twist-four analysis we have examined 
the fate of infrared renormalons. Following on from ||2^, we have demonstrated 
that infrared renormalon ambiguities appearing through the resummation of running- 
coupling effects in the twist-two coefficient functions cancel within the OPE, for both 
F2 and Fl, with another ambiguity associated with the definition of twist-four matrix 
elements of operators that mix under renormalization with twist two. While the former 
ambiguity appears due to infrared sensitivity of twist two, the latter is directly related 
to the ultraviolet renormalization properties of the twist-four operators. In fact, these 
are two manifestations of the same phenomenon, namely the arbitrariness in the sepa- 
ration of different twists. The argument that exact cancellation of the ambiguity must 
occur, is equivalent to the statement that the structure functions must be totally free of 
factorization scale dependence. 

Renormalon ambiguities pose no problem of consistency for the OPE. Full consis- 
tency to order l/Q"^ is achieved upon resumming infrared renormalons at twist two and 
including the parameterization of twist four within the same regularization prescription, 
e.g. using Borel transforms. 

The ultraviolet dominance assumption amounts to a significant simplifica- 



tion: instead of many parameters (a large number of local matrix elements), the entire 
twist-four contribution is proportional to the twist-two matrix elements, and the only 
non-perturbative parameter is the overall normalization of this contribution. Under this 
assumption, the multi-parton correlation within the target reduces to the trivial part 
which is proportional to the twist two matrix elements, whereas any more detailed cor- 
relation is neglected. From the OPE perspective this assumption seems very strong, and 
it is probably not justified for general x. On the other hand, the dominant contributions 
at large x are indeed associated with configurations similar to twist two: 

• the recoiling system corresponding to the dominant twist-four contribution to F2 at 
leading order is that of a single energetic quark (^21) . For Fl the leading-order 
recoiling system contains an energetic quark and an energetic gluon. These are the 
same configurations that appear at twist two. 
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• the initial state at twist four contains a gluon in addition to the quark or antiquark 
(in the amphtude or the complex conjugate amplitude, respectively). However, at 
large x this gluon is constrained to have a small longitudinal momentum fraction 
and, therefore, the matrix elements {p\Oa{v]y)\p) are effectively independent off. 
Whereas the detailed multi-parton correlation of twist four includes separate de- 
pendence on V and y, its dominant large-x ingredient depends only on py, just like 
the twist-two matrix elements. Consequently the number of relevant parameters 
does not increase at large x. 

Having realized that the dominant contributions at large x are associated with the 
twist-two like configurations, the assumption of ultraviolet dominance becomes natural. 
We conjecture that the simplification that occurs at large x amounts to the dominance 
of this contribution which mixes under renormalization with the leading twist. 

An additional simplification that occurs at large a:, is that certain operators, whose 
matrix element become the largest at large light-cone separations \py\ ~ iV, dominate. 
To leading order in the flavour expansion, two operators, 0^ and 03, mix with the 
leading twist. However, at large the contribution of alone dominates. Beyond the 
leading order in the flavour expansion, more operators, such as the two-gluon operators 
O5 and Oq or the four-quark operators Oi and O2 may contribute as well. 

6.2 Factorization beyond the leading twist 

Our analysis of infrared renormalons was restricted to leading order in flavour expansion 
and, in parallel, we have applied the OPE (|l^) to the Born level. It is clearly of interest 
to extend the analysis beyond this order, particularly since we know from the standard 
twist-two analysis that multi-gluon emission is important at large x. 

Higher-order corrections to the twist-two coefficient functions at large x are well 
understood |]T6|-|]19[, ^2[. They come from three different subprocesses: QCD evolution 
of the incoming quark state, q^-, hard scattering of the incoming quark off the virtual 
photon, H, and the propagation of the narrow quark jet in the final state, Jn- Due to 
different time scales involved in these three processes, they are quantum-mechanically 
incoherent, so their contribution is factorized in the moment space, 

Ff-'iN, Q') = 1 Hl{QVix') MQVNi^') g;v(/i'), (52) 

with /i^ being the factorization scale. Here, we have indicated explicitly the dependence 
of H and J on the relevant physical scales involved: the momentum transferred Q^, and 
the invariant mass of the final state = Q'^{l—x)/x ~ Q'^/N, respectively. Note that 
in the longitudinal structure function case we extracted from the hard function the factor 
so that Hl{Q'^ / fJ''^), similarly to H2{Q^/fi^) is independent of A^. The /x^ evolution 
of is induced by the emission of soft gluons with momentum k ~ Q{l—x) ^ Q/N. The 
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H dependence of the quark distribution at large follows from the renormalization 
properties of the non-local light-cone operator in the r.h.s. of (|5TD 



^-^^ = -reu.p(«s(/i')) ■ (IniV) ■ g^(/i^), (53) 

with the anomalous dimension rcusp(«s) = oisCfIt^ + C'(tts). This equation coincides 
with the DGLAP evolution equation with the matrix of anomalous dimension replaced 
by their large N asymptotic behaviour. It is important to notice that at large the 
quark distribution evolves autonomously and does not mix with the twist-two gluon 
distribution. The latter will necessary involve the emission of a soft quark into the final 
state, so the corresponding contribution is suppressed^ by a power of (1 — x) ~ 
Solving (|53D we can resum large Sudakov (double logarithmic) corrections associated 
with the infrared evolution of the quark distribution functions from the hard scale 
down to the scale Q"^ /N'^, corresponding to the transverse momenta of soft gluons in the 
final state. 

Another source of Sudakov corrections to the twist-two structure functions in (|52[) is 
the jet function J . Requiring that the l.h.s. of (|5^) should not depend on the factorization 
scale /i and taking into account that the dependence enters into the jet function only 
through as(/i^) and Q^/(A^/i^), one can obtain from (^) and ( p5[ ) the evolution equation 
on the jet function [Q. Solving this equation, one can resum large Sudakov logarithms 
to the twist-two structure functions (0). These corrections originate from an incomplete 
cancellation between virtual and real soft gluon contributions and can be attributed to 
QCD evol ution of th e outgoing quark jet carrying large energy ~ Q but small invariant 
mass ~ \JQ'^{^ — x). Let us recall |^ that the jet function depends on the nature of 
the recoiling hard partons and thus it differs in the case of the longitudinal structure 
function (a = L) with respect to the transverse one (a = 2). 

The factorization of the twist-two structure functions into the product of three factors, 
if, J and g, corresponding to three different subprocesses is a general property of QCD 
dynamics in DIS, which should be valid beyond the perturbative level. The twist four 
ingredient which mixes under renormalization with twist two has its own logarithmically- 
enhanced higher-order corrections, which include an evolution factor and a jet-function 
factor, 



Fr-\N, q\ = Hr-\QVf^') jr~\Q'/Nn') ^ gr'(/^')- (54) 

uv Q'^ 

A similar formula (with an overall supression of 1/A^) holds for F^~'^(N, Q"^) . The con- 
tribution (|54D depends, of course, on the prescription used to regularize the renormalon 
sum in J a at twist two (^21). We do not specify the regularization prescription here and 
below. The reader should keep in mind that this regularization is not associated with 
the conventional factorization (which is based on dimensional regularization and deals 

^ This can be also seen from the asymptotics of the DGLAP evolution kernel Pqg{x) (1 — x)^ as 



27 



only with the logarithmic divergence), and therefore it is not necessarily associated with 
the factorization scale fi; it can be introduced for example by the Borel transform as in 
Eq. (^). The relation between the principal value regularization of the Borel transform 
and the cutoff regularization of the corresponding momentum integrals has been inves- 



tigated in [^. The two differ just by renormalon-type power corrections of calculable 
magnitude. 

The next observation is that independence of ( [5^ ) of the factorization scale fi^ at 
the logarithmic level implies that the product of the factors H and J obeys the same 
evolution equation as the resummed twist-two coefficient function in (|52D . Thus, the 



coefficient functions in (^) differ from those in (|52|) only by the initial condition for the 
evolution 

-^iwrniwrn^' ^^"''"^^ ^^^^ 

which is a constant, = 0{N^), up to corrections suppressed by powers of as- We 
stress that nf^ is strictly independent of the ratios Q"^ j ^ and Q"^ /{Njj?). Note that n^f^ 
is not necessarily saturated by the single dressed gluon contribution, which is strictly 
associated with the light-cone operator O7. Due to contributions of the two-gluon and 
the four-quark light-cone operators, which appear beyond the large (3^ limit, /t^"'' can 
differ from of Eq. (|i2|). In particular, at this level, the power corrections to F2 and 
Fl may not be controlled by the same operators. 

In conclusion, we find that ultraviolet dominance implies that the large asymptotic 
behaviour of the structure function F2 to twist-four accuracy takes the form, 

F2(iV, Q2) = Fl^-\N, Q2) + Fl^-\N, + . . . (56) 



1 + Kf^ 



+ 



where the ellipses stand for higher-twist contributions, as well as for terms that are 
suppressed by 1/N at twist two and four. A similar factorization formula holds for 
Fl{N,Q'^). We stress that k^"^ do depend on a = 2, L. From (|55D and ( |56D it follows 
that 

4'^)A2= hm Pg^-^i^^Q") . (57) 

Here the renormalon contribution to the coefficient functions of both twist two and twist 
four is regularized. In addition, the ultraviolet divergent contribution to the twist-four 
matrix elements is regularized consistently with the renormalon regularization at twist 
two. 

Going to large x one cannot consider twist four as the only relevant power correction. 
As mentioned above, the twist six contribution is expected to scale as N"^ /Q^. This is 
reflected in the renormalon contribution |2^. The analysis of twist six within the light- 
cone expansion is signiflcantly more complicated than that of twist four. However, based 
on the experienced gained here, we assume that the same picture emerges: amongst the 
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large number of matrix elements needed to describe twist six contribution, the domi- 
nating contribution will be the one which mixes under renormalization with the leading 
twist. In particular, at Born- level for twist six, the dominant configuration to F2 will 
be the one containing a single energetic quark in the final state. Similarly, at twist 
six is dominated by a final state containing an energetic quark and an energetic gluon 
coUinear to the outgoing jet, i.e. the same configuration as the lower-twist contributions. 

Thus, "extrapolating" from the twist-four results, we expect that the dominant con- 
tributions to an?/ twist is proportional to the twist-two parton distribution g^r. This leads 
to a generalization of the large-x factorization formula ( ^61) to any order in NA'^/Q'^, 

Fl{N, Q') = ^ H^iQVlx') MQVNii') J^iNA'/Q') q^'il^'), (58) 

with 

where k„ are dimensionless non-perturbative parameters representing the leading con- 
tribution of terms with twist 2{n + 1). k^"'' was defined in (0). As argued above, the 
relative suppression of the longitudinal structure function with respect to the transverse 
one is not special to the leading twist, but is rather common to all the terms in the twist 
expansion (see e.g. (^Op). 

It thus follows from the OPE and the "ultraviolet dominance" assumption that the 
dominating power corrections in the large N limit can be all resummed into a shape 
function of a single variable J^^{NA? /Q"^) that multiplies JaiQ"^ / N fi^) . We stress that 
(^) is valid up to perturbative and non-perturbative corrections which are suppressed 
by 1/N. 

The jet function itself contains now two factors, a perturbative factor Ja{Q'^ /N ix^) 
that can be calculated simply using the twist-two operators, and a non-perturbative 
factor J^^{NK^/Q'^) which sums up power corrections to all orders. The separation 
between them by regularizing the renormalon sum in JaiQ"^ / N fi"^) is, of course, arbitrary. 
However, the product is regularization independent: both factors together describe the 
hadronization of the recoiling jet. Since the hadronization process takes place through 
the exchange of soft gluons with the remnants of the target, the parameters k^"-*, and 
thus the function J^^ {N A'^ / Q"^) , do depend]^ on the target. 

Our analysis in this paper was restricted to the parity-conserving spin-averaged struc- 
ture functions F2 and Fl. If the ultraviolet dominance conjecture indeed holds, it is nat- 
ural to expect that large- iV factorization formulae such as ( ^8]) apply to a wider class of 

^"^In ref. the factor Ja is defined in a target independent way. On the other hand the power correc- 
tions encoded by J^^ contain information about the matrix elements of higher twist operators, through 
the target dependent constants k^^^ . Thus we cannot interpret J^^ as merely the non-perturbative cor- 
rections to Ja, but rather as a more general, target dependent factor which encodes the entire power 
corrections to the moments of the structure functions. 
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structure functions. In particular, we expect it to apply to the parity-violating structure 
function F3 as well as the spin structure function gi, which both share ||2^ the same 



leading twist 1/Q^ renormalon ambiguity with Fi = ^(-^2 — Fl). Note, however, that 
in order to identify the relevant twist-four operators and trace down the cancellation of 
renormalon ambiguities one must first work out the operator expressions for the twist 



four contribution, similarly to the ones of |]rT|, which are not yet available for these 
structure functions. 

Let us try to understand the physical interpretation of the obtained expression for 
the structure functions (pSf). It is convenient to introduce a new non-perturbative dis- 
tribution function jaifn"^) of a single argument [| 
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Jo 



(60) 



By the definition, this function is an inclusive distribution describing the fragmentation 
of a partonic system with the invariant mass into the hadronic final state, jaifn"^) 
does not depend on the momentum transfer and, it is, in this sense, a universal 
distribution function. Nevertheless, it does depend on the target as well as on the 
renormalon regularization prescription (see below). Inserting (^Dp into ( |55D and going 
over from the moment space to the momentum-fraction representation]^ we obtain the 
following expression for the structure function with the power corrections included. 



J X 



2,L), (61) 



where we assumed that jaim'^) vanishes fast for ^ Q^(l — x). Here g^^^^jg ^qj- 

the twist-two structure function in which infrared renormalon in the resummed coefficient 
functions have been regularized. 

We conclude that the effect of the power corrections to all orders in {N/Q"^) is to 
smear the twist-two structure functions close to the end-point, x ~ 1. Whereas both 
and ja depend on the renormalon regularization prescription, this dependence 
cancels out in their convolution in the r.h.s. of Eq. (pT]). 

Let us consider the general properties of the distribution function ja(m^). The OPE 
based relation ( |59D leads to a similar expansion of the function for ja(m^). 



+ ■ 



(62) 



This expansion is singular as it goes over the derivatives of the delta function. Its 
substitution into ( pl| ) leads to the following twist expansion at large x 



Fa{x,Q' 



-itw— 2 J 



(a) 



g2(i-x) 



:i-x)^Fr'-\x,Q') 



+ K,2 



ii-xr-f,F:^''i 



+ 



(63) 



^^At the accuracy considered, i.e. when terms suppressed by 1/iV are neglected ^ ~ e ^\ so 
the Mellin transform is equivalent to a Laplace transform. 
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where higher-twist corrections are proportional to the derivatives of twist two. The 
twist-four correction agrees with the expression obtained in Ref. p9|.P^ 

In order to reconstruct the shape of the distribution ja(m^), one has to resum the 
whole series in N/Q"^ in (^). Obviously, this will require a full knowledge of non- 
perturbative effects, which is far beyond reach. Nevertheless, based on a general prop- 
erties of non-perturbative distributions [Q, we expect that ja(^^) should vanish for 

< 0, increase as a power at small m^, approach its maximal value at ~ and, 
then, rapidly decrease at large ~ — x). 

Moreover, there is a way to probe the structure of the non-perturbative distribution 
jaifn^) by Dressed Gluon Exponentiation [Q, deducing an ansatz for Ja{m'^). Rely- 
ing on the perturbative treatment of soft and collinear radiation one can conclude that 
renormalon related power-corrections exponentiate as well. This can be systematically 
established by studying the large-order behaviour of the Sudakov exponent: it turns out 
that the coefficients of sub-leading Sudakov logarithms are enhanced factorially due to 



infrared renormalons The enhancement of sub-leading logs implies that Sudakov 
resummation with a fixed logarithmic accuracy is insufficient in the limit under consid- 
eration. Instead, one is obliged to resum all the logarithms that are factorially enhanced 
in Ja- As usual, when is calculated perturbatively to power accuracy, a renormalon 
ambiguity shows up. This ambiguity is compensated at the non-perturbative level by 
J^^. The fact that the renormalon ambiguity in the Sudakov exponent is additive is 



fully consistent with the factorized form of (|58D we obtained above. 



At the level of a single dressed gluon, the exponent is equal to the logarithmically- 
enhanced part in the single dressed gluon result. It therefore has two Borel singularities, 
at M = 1, 2, suggesting the following ansatz for the non-perturbative jet function P5[ , 



J„ (A^A IQ ) = exp - -4 ^^^1 ^ (64) 

where Wi'^'* = — k^"^ and = —i^^2^ + \ ('^1"^) determine to the center and the width 
of the distribution ja{fn^)i respectively. Clearly, c^i 2 depend on the regularization of the 
renormalon sum in the resummed coefficient function J a- It is interesting to note that 
the same structure of the exponent, and thus the same ansatz for power corrections in 
the large-x region, appears also in the case of the transverse and the total fragmentation 
functions in e^e~ annihilation. 

Eq. (|6^) should be regarded as a minimal ansatz for the parametrization of power 
corrections at large x, since further renormalon ambiguities in will probably appear 
once the analysis of the Sudakov exponent is generalized to include more than a single 



dressed gluon. In this case further power corrections should be included in the exponent 
of Eq. (|6^). We must stress, that the suppression of such multi-gluon contributions to the 
non-perturbative exponent is not at all obvious, and it remains to be checked]^. On the 

^^We are grateful to J. W. Qiu for drawing to our attention to this paper. 

^•^Some phenomenological evidence of such suppression was found in in the case of the thrust 
distribution in e+e~ annihilation, barring hadron mass corrections. 
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other hand, a non-perturbative distribution which is more involved than (0) can only 
be consistently used provided that renormalon resummation in at the corresponding 
level is performed. 
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